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Abstract 

We study the LHC search prospects for a model in which the neutrinos obtain Dirac masses from 
couplings to a second Higgs doublet with tiny vacuum expectation value. The model contains a 
charged Higgs boson that decays to £v with branching fractions controlled by the neutrino masses 
and mixing angles as measured in neutrino oscillation experiments. The most promising signal 
is electroweak production of H + H~ pairs with decays to ££'p™ lss , where ££' = e + e~, u + u~, and 
e ± /i T . We find that a cut on the kinematic variable Mt2 eliminates most of the tt and VK-pair 
background. Depending on the neutrino mass spectrum and mixing angles, a 100 (300) GeV 
charged Higgs could be discovered at the LHC with as little as 8 (24) fb _1 of integrated luminosity 
at 14 TeV pp center-of-mass energy. 
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I. INTRODUCTION 



The Standard Model (SM) accounts for almost all experimental high energy physics 
data; however, the observation of neutrino oscillations requires that the SM be extended to 
include nonzero neutrino masses. While there are many ways to expand the SM to account 
for neutrino oscillations, we attempt to do so with the following goals. First, the neutrino 
mass scale is significantly lower than the mass scales of the other fermions, so we would like 
the model to account for this without the addition of many tiny parameters. Second, lepton 
number violation has not yet been observed, so we would like the model to give rise to Dirac 
neutrino masses, with Majorana masses forbidden. Third, we would like the model to be 
testable at the CERN Large Hadron Collider (LHC). 

Most neutrino mass models give rise to Majorana masses for the SM neutrinos, with many 
predicting TeV-scale new physics accessible at the LHC. In contrast, only a few models for 
Dirac neutrinos have been proposed. These typically involve a second Higgs doublet with 
very small vacuum expectation value (vev) that couples only to the left-handed lepton 
doublets and the right-handed neutrinos, resulting in neutrino masses of the same order as 
the very small vev. The original SM-like Higgs doublet couples to all of the quarks and 
charged leptons in the usual way. Such a Yukawa coupling structure can be obtained by 
imposing a global Z 2 symmetry, as proposed in the models of Refs. jl, 2|; however, this does 
not by itself forbid neutrino Majorana mass terms, which must instead be eliminated by 
imposing an additional lepton number symmetry. The required Yukawa coupling structure 
can also be obtained by imposing a global U(l) symmetry; this idea was first proposed in 
Ref. [3 1 as a way of ensuring the (then-assumed) masslessness of the neutrinos in the presence 
of right-handed neutrino states, and has the virtue of forbidding Majorana mass terms by 
itself. 

In order to generate neutrino masses, the global symmetry used to ensure the desired 
Yukawa structure has to be broken. Spontaneous breaking leads to a very light scalar which 
can cause problems with standard big-bang nucleosynthesis Q , as well as having significant 
effects on the phenomenology of the new Higgs particles 0|. By instead breaking a global 
U(l) symmetry explicitly, the model proposed by us in Ref. [5] generates Dirac neutrino 
masses while avoiding very light scalarsjj A supersymmetric version of this model was 



1 A similar mechanism was used to explain the top-bottom quark mass hierarchy in Ref. 
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studied in Ref. which found spectacular multi-lepton signals from cascade decays of the 
supersymmetric partners of the new Higgs bosons and right-handed neutrinos at the LHC. 

In this paper we study the LHC detection prospects of the non-supersymmetric model 
of Ref. jsj. This model expands the SM by adding a second Higgs doublet $2 with the 
same electroweak quantum numbers as the SM Higgs doublet $1, as well as adding three 
gauge-singlet right-handed Weyl spinors Vr. which will become the right-handed components 
of the three Dirac neutrinos. The model imposes a global U(l) symmetry under which the 
second Higgs doublet and the right-handed neutrinos have charge +1, while all the SM fields 
have charge zero. This allows Yukawa couplings of the second Higgs doublet only to the 
right-handed neutrinos and the SM lepton doublet, and forbids Majorana masses for the 
right-handed neutrinos. It also tightly constrains the form of the Higgs potential. Breaking 
the U(l) symmetry explicitly using a term m\ 2 ^\^2 + h.c. in the Higgs potential yields 
a vev v 2 for the second Higgs doublet and consequently gives the neutrinos Dirac masses 
proportional to v 2 - By requiring that v 2 ~ 0(eV), the Dirac neutrino masses are made 
suitably small without requiring tiny Yukawa couplings. 

The characteristic feature of the model is that the couplings of the charged scalar pair if ± 
and two neutral scalars H° and A from the second Higgs doublet to leptons and neutrinos 
are controlled by the neutrino masses and mixing angles. In this paper we take advantage of 
the distinctive decay of the charged Higgs boson H + into charged leptons and neutrinos. We 
focus on electroweak pair production of H + H~ at the LHC followed by decays to ££'p™ lss , 
where ££' can be any combination of opposite-sign e, /x, and r leptons and p™ lss denotes 
missing transverse momentum (carried away by the neutrinos). Because r leptons are more 
difficult to reconstruct experimentally, we concentrate on the final states with ££' = e + e~, 
and e ± fi T . The major backgrounds are diboson production (W + W~, ZZ, and Zj) 
and top quark pair production with both tops decaying leptonically. 

To determine whether the H + H~ signal will be detectable at the LHC, we generated 



signal and background events using MadGraph/MadEvent version 4 js| assuming 14 TeV pp 
center-of-mass energy. We present results both at parton level, and after hadronization with 



PYTHIA [9] and fast detector simulation with PGS 10|. With appropriate cuts, we find that 
a 5cr discovery can be achieved with luminosity in the range 8-75 fb _1 for M H + = 100 GeV, 
depending on the neutrino mixing parameters. For Mh+ = 300 GeV a 5er discovery can be 
made with luminosity in the range 24-460 fb _1 . The higher luminosity requirements occur 



when the neutrino parameters are such that H + decays mostly to ru, leading to final states 
not considered in our analysis. We find that the kinematic variable Mt2 is very effective at 
separating the signal from the ti and WW backgrounds for charged Higgs masses above the 
W mass, and also provides sensitivity to the charged Higgs massif 

This paper is organized as follows. In the next section we review the model and present 
the charged Higgs decay branching ratios. In Sec. IHII we describe the signal and background 
processes, our event generation procedure and selection cuts, and the resulting signal signif- 
icance. In Sec. IIVI we summarize our conclusions. 

II. THE MODEL 

As outlined in the introduction, we start with the field content of the SM and add to it a 
new scalar SU(2)l doublet $2 (the SM Higgs is denoted $1) and three right-handed gauge 
singlets Vr. (these are the right-handed neutrinos). We impose a U(l) symmetry under 
which $2 an d the three have charge +1 and all the other fields are uncharged, which 
leads to the Yukawa coupling structure j^] 



Here $j = icr 2 $* is the conjugate Higgs doublet and y{- are the 3x3 Yukawa matrices for 
fermion species /. 

The Higgs doublets can be written explicitly as 



where v\ will be generated by the usual spontaneous symmetry breaking mechanism of the 

SM and v-i will be generated by the explicit breaking of the global U(l), described below. 

Inserting these expressions for $j into Eq. (JTJ, we obtain the fermion masses and couplings 

to scalars. In particular, the fourth term in Eq. ([TJ gives rise to the neutrino mass matrix 

2 While we have not made a detailed study of charged Higgs detection prospects at 7 TeV pp centre-of- 
mass energy, we note that the cross section for the most dangerous WW background is about 2.5 times 
smaller at 7 TeV. However, the signal cross section is also about 2.5 (4.5) times smaller at this energy for 
M H + = 100 (300) GeV. Furthermore, the LHC is anticipated to collect only about 1 fb _1 of integrated 
luminosity at 7 TeV. We thus expect detection or even exclusion of the process considered here to be 
unfeasible in the current 7 TeV LHC run. 



yfjduMQL, ~ V^urMQl, - vlfn^XL^ - y^v R ^\L Lj + h.c. 



(1) 
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and interactions: 



r _ Uij v 2 _ Vij ,0 r- -Hii iO,i- . v ± + - n i l_ 



(3) 



After diagonalizing the mass matrix in the first term, the neutrino mass eigenvalues are 
given by m Vi = y^v 2 /V2, where y\ are the eigenvalues of y^. In this way, the small masses 
of the three neutrinos can be traced to the small value of v 2 . 

We obtain the vevs of the scalar doublets from the Higgs potential as follows. The most 
general gauge- invariant scalar potential for two Higgs doublets is (see, e.g., Ref. 11]). 



V = mj^l^i + m 2 22 §\§ 2 - mi 2 $|$2 + h.c. 

+^Ai ($1$ x) 2 + ^A 2 ($J$ a ) 2 + A 3 ($ ($J$ 2 ) + A 4 ($I$ 2 ) 



+ \ ^A 5 ($I$ 2 ) 2 + A 6 $I$i + A 7 $ 2 $ 2 $l$ 2 +h.c. 



(4) 



Imposing the global U(l) symmetry eliminates m\ 2l A5, Ag, and A7. The global U(l) sym- 
metry is broken explicitly by reintroducing a small value for m\ 2 . This leaves the Higgs 
potential j^f] 



1 



mi 2 ${$2 + h.c. 



+ ^Ax (V^x)' + ^A 2 ($ 2 $ 2 ) 2 + A 3 (V^i) ($J$ 2 ) + A 4 ($J$ 2 ) 



(5) 



Stability of the potential at large field values requires A 1; A 2 > 0, A 3 > — \A1A2, and A4 > 
— vAiA~2 — A3. We want v\ to arise through the usual spontaneous symmetry breaking 
mechanism, which is achieved when vn? xl < 0. We do not want the global U(l) to also be 
broken spontaneously, as that will create a very light pseudo-Nambu-Goldstone boson, which 
is incompatible with standard big-bang nucleosynthesis; thus we require that the curvature of 
the potential in the v 2 direction at zero <3> 2 field value be positive, i.e., m\ 2 + {\^+\4)v 2 /2 > 0. 

To find the values of the vevs in terms of the parameters of the Higgs potential, we apply 
the minimization conditions, 



dV 



m 2 xl vi - m\ 2 v 2 + + o^ 3 + ^) VlV 







Note that using a Z2 symmetry instead of the global U(l) would allow a nonzero A5 term. 
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dV 



m\ 2 v 2 - m\ 2 v x + -A 2 i| + ^(A 3 + A 4 )^ 2 = 0. (6) 



Since we will require m\ 2 C fj, we can ignore m\ 2 and v 2 when finding the value of v\ . This 
yields 

—2m 2 

v\ = -^li. (7) 



For t> 2 , we need to consider mf 2 , although again we may ignore higher order terms in m\ 2 jv\\ 
this yields 

m\ 2 v\ 

V2 = m| 2 + i(A 3 + A 4 K- (8) 

We will choose parameters so that v i ~ 246 GeV and t> 2 ~ eV. This requires m\ 2 ~ (MeV) 2 . 
We note that because m\ 2 is the only source of breaking of the global U(l) symmetry, its 
size is technically natural; i.e., radiative corrections to m\ 2 are proportional to m\ 2 itself 
and are only logarithmically sensitive to the high-scale cut-off 



The mass eigenstates of the charged and CP-odd neutral scalars are given by 

G + = <f)f sin p + <\>% cos ~ (f>i 
H + = <f)f cos ft — (j) 2 sin ft ~ — cj) 2 
G° = # sin/3 + 0^ cos/3~0^ 

A Q = ^cos/3-0^sin^~ (9) 

where we define tan/3 = t>i/t> 2 ~ 10 11 . G + and G° are the Goldstone bosons, which do not 
appear as physical particles in the unitarity gauge. H + and A are the physical charged and 
CP-odd neutral Higgs states and are almost entirely contained in $ 2 . Neglecting contribu- 
tions of order m\ 2 and v 2 , the masses of H + and A are 

M 2 H+ = m\ 2 + ^X 3 vl 
M\ = m\ 2 + i(A 3 + A 4 H 2 = M 2 H+ + ^4^. (10) 

The mass matrix for the CP-even neutral states is almost diagonal, yielding only very tiny 
mixing of order v 2 /v\. Ignoring the mixing, the eigenstates are h° ~ 4>Y (SM-like) and 
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Parameter 


Value 


sin ^ f)i n 

Oil! l/j_2 

sin 2 9 23 
sin 2 


314fl + nit N ) 

U - Ji ^l i -0.15J 

0.44(li^) 
0.9± 2 ; 3 x 10~ 2 


Am 2 = m 2 2 — m 2 x 


7.92(1 ±0.09) x 10~ 5 eV 2 
±2.4(l^- 2 i) x 10" 3 eV 2 



TABLE I. Current values of the neutrino mixing parameters and mass-squared differences, from 
the global fit to neutrino oscillation data performed in Ref. 12] . Uncertainties quoted are the 2a 
ranges. Note that the constraint on sin #13 is only an upper bound, and that the sign of AM 2 is 
not yet known. 



~ (f>2 ,r , with masses 



Ml 



M: 



H 



m 22 + ^( A 3 + A 4 )^ 2 ■ 



Ml 



(11) 



After diagonalizing the neutrino mass matrix, Eq. (]3]) yields the following couplings to 
the new physical Higgs states: 



^Yuk ^ 



m 



H*m + i^A'Pa^ - y^lUlH+u^et + h.c], 
v 2 v 2 v 2 



(12) 



where Un is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, defined according to 
Vi = Ua^i, where V(_ are the neutrino flavor eigenstates. 

The PMNS matrix can be parameterized in terms of three mixing angles dij (with ij = 12, 
23, and 13) and a phase 5 according to (see, e.g., Ref. jl2|). 



( 



\ 



C12C13 



S12C13 



s\ 3 e~ 



\ 



-S12C23 - C12S23S13C" 5 C12C23 - S12S23S13C" 5 S23C13 

S12S23 - ci 2 c 23 s 13 e tS -c 12 s 23 - s 12 c 23 s 13 e tS c 23 c 13 



J 



(13) 



where c,- 



cos 6ij and s^- 



sin^j. The 2a experimentally-allowed ranges for the three 
mixing angles and the neutrino mass-squared differences are given in Table [U The phase 
5 and the mass of the lightest neutrino are undetermined, although tritium beta decay 
experiments set an upper limit on the neutrino masses of about 2 eV [l3l ]. 

Since the decays of H° and A to two neutrinos will be invisible to a collider detector, the 



decay of most interest is H + — > £ + v. The charged Higgs can decay into all nine combinations 
of £ii>j] summing over neutrino mass eigenstates, the partial width to a particular charged 
lepton £ is [5| 

r -+ £+u) = M *+<rfh (14) 

where we define the expectation value of the neutrino mass-squared in a flavor eigenstate 
by (l4| 

i 

In what follows we work under the assumption that Mga j\o > Mh+, i.e., A4 > 0, so that 
the decays H + — > W + H°, W + A° will be kinematically forbidden. The branching ratios of 
the charged Higgs are then completely determined by the neutrino masses and mixing: 

where we used the unitarity of the PMNS matrix to simplify the denominator. 

The sign of the larger neutrino mass splitting AM 2 is unknown (see Table [TJ). The 
situation in which AM 2 is positive, so that z/3 is the heaviest neutrino, is called the normal 
neutrino mass hierarchy, while the situation in which AM 2 is negative, so that v\ and v 2 are 
heavier, is called the inverted hierarchy. We compute the charged Higgs branching fractions 
as a function of the lightest neutrino mass for both hierarchies, scanning over the 2a allowed 
ranges of the neutrino parameters as given in Table [TJ Results are shown in Fig. [Tin The 
large spread in the branching ratios to \iv and rv for lightest neutrino masses below about 
0.06 eV is due to the current experimental uncertainty in sin 2 #23, which controls the relative 
amount of and v T in the isolated mass eigenstate z/ 3 . 

1 — I 

Limits on the model parameters were discussed in Ref. [5|]. The most significant for our 
purposes is from searches for leptons plus missing energy at the CERN Large Electron- 
Positron Collider, which put a lower bound on the charged Higgs mass of 65-85 GeV, 
depending on the mass of the lightest neutrino. Big-bang nucleosynthesis also puts an 



4 We disagree with the charged Higgs branching fractions to leptons presented in Ref. [4j for the Zi model 
of Ref. [2| ; these decays should have the same relative branching fractions as in our model. 
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FIG. 1. Charged Higgs decay branching fractions to ev, p,v, and tv as a function of the lightest 
neutrino mass. 



upper bound on the neutrino Yukawa couplings of 



V-2 



~ 30 



M H+ 
100 GeV 



\U, 



it 



(17) 



III. SIGNAL AND BACKGROUND AT THE LHC 



In most other two-Higgs-doublet models, the charged Higgs decay rate to a particu 



ar 



charged lepton is proportional to the square of the charged lepton mass (see, e.g., Ref. [11]). 
Such a charged Higgs therefore decays predominantly to tv, with decays to \iv, ev below 1%. 
In our neutrino-mass model, however, the charged Higgs decay rate to a particular charged 
lepton is instead proportional to the square of the mass of the corresponding neutrino flavor 
eigenstate. As a result, the branching fraction to ev and/or [iv will always be sizable. 
In particular, in the normal hierarchy BR(iJ + — > \w) ~ 1/2, in the inverted hierarchy 
BK(H + — > ev) ~ 1/2 and BR(iJ + —>• \xv) ~ 1/4, and for a degenerate neutrino spectrum 
BR(iJ + — > ev) ~ BR(if + — > fiv) ~ 1/3, as shown in Fig. [U Considering the high detection 
efficiency and lower fake rates of e and /i compared to r, we study H + H~ pair production 
at the LHC mediated by a photon or Z, followed by decays to e or fi with missing transverse 
momentum. We consider two scenarios, = 100 and 300 GeV, and present results as a 

function of the lightest neutrino mass for both the normal and inverted hierarchy. 
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Rranrhinp* frartinn 

1^1 CAillV^lllll!-. J--L CAil^/ U J.v_yJ.J. 


W + £+z/ £ (£ = e or /x) 


0.1068 


Z -> (£ = e or /i) 


0.0336 


Z vv (all 3 neutrinos) 


0.2000 




1.0000 



TABLE II. Default SM branching fractions used in MadGraph/MadEvent @]. 

The process of interest is pp —¥ H + H~ — > ii'&tisp, with = e + e~, and e ± /i T . 

The relevant backgrounds are — >■ VV — )■ ^'Pv with VV = W + W~, ZZ, or Z'-f and the 
neutrinos of any type, and pp — > ti — > ££'&iVj»rbb. In spite of the presence of the extra b 
jets that can be vetoed, the ti process is important because of its exceptionally high cross 
section at the LHC. 



A. Event generation 

We simulated the signal and background processes with the parton-level Monte Carlo 
MadGraph/MadEvent version 4 [8J]. We present both a parton-level analysis and an analysis 
including showering, hadronization, and a fast detector simulation using a PYTHIA-PGS 
package designed to be used with MadEvent. PYTHIA (version 6.4.20) [9(] generates initial- 
and final-state radiation and hadronizes the final-state quarks and gluons, while PGS (Pretty 
Good Simulation of High Energy Collisions, version 4) 10( is a basic detector simulator — we 
used the default settings for ATLAS. For the signal process we generated 10,000 unweighted 
events in each of the e + e~, /i + p~, and /i + e - final states. For both the VV and ti backgrounds 
we generated 100,000 unweighted events in each of the three leptonic final states. We 
incorporated the p~e + final state by doubling the fi + e~ cross sections. For the backgrounds 
we used the default SM branching fractions from MadGraph/MadEvent, given in Table HT1 

Although MadGraph/MadEvent is a tree- level event generator, we partially incorporated 
next-to-leading order (NLO) QCD corrections. We did this for two reasons. First, QCD 
corrections have a significant effect on the signal and background (especially tt) cross sec- 
tions, as well as significantly reducing the QCD scale uncertainty, so that using NLO cross 
sections lets us obtain more reliable results. Second, for the M H + = 100 GeV simulation 
we will apply a jet veto to reduce the ti background, which will also affect the signal and 
VV background once initial-state radiation is included. While this could be simulated by 
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running the no-jet events through PYTHIA, a parton-level simulation of the H + H~j and 
VVj processes provides a more accurate description of jet kinematics. Because these one-jet 
processes make up part of the NLO QCD cross section for the corresponding no-jet processes, 
we must incorporate the NLO cross sections for consistency, as follows. 

In the absence of a full NLO Monte Carlo, NLO QCD corrections are usually incorporated 
by multiplying the leading-order (LO) cross section — and the cross section corresponding to 
each simulated event both before and after cuts — by a /c-factor equal to the ratio of the NLO 
cross section to the tree-level cross section. In our case, however, our jet veto will affect LO 
events (which have no jet) and NLO events (which can have a final-state jet) differently. 
We deal with this by simulating pp — > H + H~j and pp — > VVj with the same decay final 
states as considered in the no-jet processes. For simplicity we generate the same number of 
events with an additional jet at the parton level as were generated for the no-jet processes. 
Because the ti background already contains two jets at leading order, we do not separately 
generate events with additional jets for this background. To avoid the collinear and infrared 
singularities, we apply a minimum pt cut of 10 GeV on the jet at the event-generation level. 

The square of the NLO matrix element can be expressed up to order a s as 

IMInlo = \Mjjo + -M lloop | 2 + |M|; jet . (18) 
We used MadGraph/MadEvent to calculate the cross sections corresponding to A^lo and 



M ijet- We computed the N 
FORTRAN code PR0SPIN0 



O cross-section for pp — > H + H at the 



15 



;HC using the public 



16] with CTEQ6 parton densities 17|, with the renormal- 
ization and factorization scales set equal to M#+. We took the NLO cross sections for the 
SM W + W~ and ZZ background processes from Ref. [la]. This paper quotes results using 
both the MRS98 and CTEQ5 parton densities, with results differing by ~6%; since we use 
CTEQ6 for the tree-level MadGraph/MadEvent calculation, we take the results using the 
CTEQ5 parton densities for consistency. For events with e ± fi T in the final state, only the 
cross section for W + W~ is relevant; for events with fi + fi~ or e + e~ in the final state, both 
the W + W~ and ZZ processes contribute and we add the cross sections at both LO and 



NLO. We took the ti cross section from Ref. [19|], which includes both NLO and next-to- 
leading logarithmic corrections. The remaining scale uncertainty is about ±5% when the 
factorization and renormalization scales are varied between m t /2 and 2m t . The relevant 
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Process 



Cross section Source 



pp H+H- (M H+ = 100 GeV) 295 fb PR0SPIN0 [15, 16] 
pp -> iJ+iJ- (M^+ = 300 GeV) 5.32 fb PR0SPIN0 [15, 16] 
W + W~ 127.8 pb Ref. [18] 
pp -»■ 17.2 pb Ref. [18] 
-> tt 833 pb Ref. [19] 



TABLE III. NLO cross sections for signal and background processes (before decays) at the LHC 
(14 TeV). The tt cross section also includes a resummation of next-to-leading logarithmic correc- 
tions. 



cross sections are given in Table IHIl 

We find that with our generator-level jet pt cut on <7ij et , cxnlo < °lo + ^ljet, so the 
one-loop matrix element must interfere destructively with the LO matrix element. Thus the 
generated cross section from the LO process must be scaled down in order to incorporate the 
effects of the one-loop correction. For the parton-level simulation, the relevant scale factor 
is determined by solving for k in the equation, 

onlo = ka-Lo + 0i jet, (19) 

before cuts are applied, and then using this equation with the same value of k to calculate 
the surviving ctnlo after the cuts are applied to the LO and one-jet MadGraph/MadEvent 
simulated results. 

For the PYTHIA-PGS simulation, the simulated events have extra jets produced by 
PYTHIA and "measured" jet px smeared by PGS. To avoid double-counting, we use the 
following equation with two constants: 

a NLO = ma^ + na ( [ f ct , (20) 

where <t£o and erf^t are the cross sections identified by PGS as having no jets and at least 
one jet, respectively, with px > 10 GeV, out of the combined LO and one-jet generated 
samples. The constants m and n are determined by maf^ = ka^o and ncr^ = crijet using 
k from Eq. (|T9|) . Equation ( 120]) with the same values of m and n is then used after cuts to 
calculate the surviving <7nlo- 

12 



Cut name Explanation 

Basic cuts Present are a lepton and antilepton, each with p T > 20 GeV, and missing 
transverse momentum p™ lss > 30 GeV. For the parton level results, we 
also apply lepton acceptance cuts of \q e \ < 3.0 and |r/^| < 2.4. 

Z pole veto To eliminate events that include Z — > £ + £~ , we veto events in which the 
invariant mass of e + e~ or fi + fi~ is between 80 and 100 GeV (not applied 
to the e ± fi T p T llss final state). 

Jet veto Designed to reduce tt background, any event with a jet with p^ > 
30 GeV was rejected. For the parton level results, this veto is only 
applied when |?7j C t| < 5.0. (Applied only for M H + = 100 GeV.) 

M T2 cut To reduce the W + W~ and tt backgrounds, we make use of the larger 
mass of H + compared to the intermediate W bosons in both backgrounds 
by cutting on M T2 (defined in Eq. 021])). For M H + = 100 GeV we 
require M w < M T2 < 100 GeV and for M H + = 300 GeV we require 
150 GeV < M T2 < 300 GeV. 



TABLE IV. Summary of cuts. 



B. Cuts 



We apply four cuts to reduce the background, summarized in Table [TVl The first cut 
checks for the presence of two opposite-sign leptons each with px > 20 GeV and missing 
transverse momentum of at least 30 GeV. For the parton-level simulation, we also apply 
acceptance cuts on the pseudorapidity of both leptons, \r)\ < 3.0 for electrons and \r)\ < 2.4 
for muons. Second, for the e + e~ and final states we veto events for which the dilepton 

invariant mass falls between 80 and 100 GeV, in order to eliminate background from Z{— > 
££) + p™ lss . This will also eliminate the majority of any background from Z + jets with 
fake p™ lss , which we did not simulate. The third cut vetoes events containing a jet with 
Pt > 30 GeV; for the parton-level simulation, we require that this jet falls in the rapidity 
range \rj\ < 5.0. This eliminates more than 97% of the ti backgound, but also reduces the 
signal by about a factor of two. We find that this cut is useful for M#+ = 100 GeV. For 
M H + = 300 GeV the signal cross section is considerably smaller and the signal events will 
be better separated from background in our final cut variable, so that we obtain better 
sensitivity without the jet veto. 
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The final cut is on the variable Mt2, defined as 



20] 



_ mm 

ML = 

1 * missfl) . miss(2) m \ 

q T + q T = 



, 2 ( t(l) missfl) \ 2 ( €(2) miss(2) 

max < m T I p 7 : , q T I , m T I p 7 : , g T 



(21) 



where is the square of the transverse mass (ignoring the charged lepton and neutrino 
masses) , 

m% (p T , q^ iss ) = 2 (\p T \ |g™ ss | - p* • <f™ ss ) . (22) 

In other words, Mt2 is determined by making a guess for the transverse momenta of the two 
neutrinos (constrained by the measured total missing transverse momentum) and computing 
the transverse masses of the two Ev systems; the guess is then varied until the larger of the 
two reconstructed transverse masses is minimized. 

For equal-mass intermediate particles each decaying to £u, the Mt2 distribution has an 
upper endpoint at the mass of the intermediate particle. Thus by cutting out events with 
Mt2 < Mw, all the W + W~ background should be eliminated (the endpoint is in fact 
smeared out by the finite W width and momentum resolution of the detector). Since the 
leptons and missing transverse momentum in the ti process also come from decays of on-shell 
W + W~, this background should be eliminated as well. There is also a small contribution to 
the VV background from nonresonant processes that can have Mt2 > Mw- Since all signal 
events will have Mt2 < Mjj+ , we also cut out events with Mt2 > Mh+ in an effort to reduce 
the background from these nonresonant VV processes. For Mh+ = 300 GeV, we find that 
raising the minimum cut on My 2 to 150 GeV reduces the tail of the nonresonant VV events 
without reducing the signal too much. 

In Fig. [2] we show the M^ 2 distributions for signal and background processes in the 
e + e~p miss channel for M#+ = 100 and 300 GeV after the other cuts have been applied, for 
the PYTHIA-PGS simulation. Note that the tt background distribution has a maximum 
Mt2 value a little above the W mass, so that it can be eliminated with a high enough cut 
on Mt2, as we do for the case of M#+ = 300 GeV. (The higher M^2 endpoint for tt in the 
right-hand plot in Fig. |2] is due to the absence of the jet veto, resulting in much higher tt 
statistics.) The VV background also falls off dramatically around M T2 ~ M w ; however, due 
to nonresonant diagrams without on-shell intermediate W pairs, this background extends 
to much higher values of Mt2- With our simulation statistics, a single tt event corresponds 
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Mj 2 [GeV] Mj 2 [GeV] 

FIG. 2. Mt2 distributions after other cuts have been applied for the e + e~p™ lss final state, with 
M H + = 100 GeV (left, with jet veto) and 300 GeV (right, no jet veto). For the signal we take 
BR(i7 + —J- e + u) = 1/3, which occurs for a degenerate neutrino spectrum. The Mt2 cut window is 
shown by the vertical lines. 

to a cross section of about 0.1 fb, while a single VV event corresponds to a cross section of 
about 0.01 fb. 

C. Results 

The efficiency of each cut on ctnlo for the e + e~p™ lss final state is displayed in Tables PVT. 
ED and IVHI Cut efficiencies for /i + /i"p™ lss are displayed in Tables IVIH[ HXl and [Xj, and for 
e ± /i =F p™ lss in Tables IXI] IXIII and lXIIII We give efficiencies for both the parton-level simulation 
and the simulation including showering, hadronization, and fast detector simulation using the 
PYTHIA-PGS package. All results incorporate NLO corrections as described in Sec. IIII A[ 

Consider for example the PYTHIA-PGS results in the e + e~p^ lss final state, and assume 
a degenerate neutrino spectrum so that BR,(H + — > e + v) = 1/3. In this case, for M#+ = 
100 GeV, the cuts reduce the charged Higgs signal cross section in this channel from 32.8 fb 
to 1.34 fb, while reducing the VV background from 1570 fb to 3.01 fb and the tt background 
from 9500 fb to 2.57 fb. The ratio of signal to background cross sections (S/B) is then 
0.24. For Mh+ = 300 GeV, S/B is comparable. These are displayed for all channels for a 
degenerate neutrino spectrum in Table IXTVl In all cases S/B is at least 0.22, comfortably 
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= 100 GeV 


M H+ 


= oUU (jreV 


Cuts Parton PYTHIA/PGS Parton 




Basic cuts 0.62644 


0.47860 


0.92961 


72105 


Z pole veto 0.90754 


0.90165 


0.97732 


0.97724 


Jet veto 0.68433 


0.60717 






M T2 cut 0.17075 


0.15542 


0.47317 


0.45945 


Cumulative 0.06643 


0.04072 


0.42988 


0.32375 



TABLE V. Cut efficiencies for the signal process pp — > e + e~p^ lss via H + H~ . The efficiency of 
each cut is defined as the cross section that passed the cut divided by the cross section that passed 
the previous cut. The cumulative efficiency is the cross section that passed all the cuts divided 
by the original cross section. The jet veto is not applied for M H + = 300 GeV. The Mt 2 cut is 
M w < M T2 < 100 GeV for M H + = 100 GeV, and 150 GeV < M T2 < 300 GeV for M H+ = 300 GeV. 





VV Background 


tt Background 


Cuts 


Parton 


PYTHIA/PGS Parton 


PYTHIA/PGS 


Basic cuts 


0.42708 


0.33912 


0.58407 


0.40612 


Z pole veto 


0.74727 


0.73255 


0.86236 


0.85501 


Jet veto 


0.63306 


0.67299 


0.01318 


0.02856 


M w + < M T2 


< 100 GeV 0.01401 


0.01147 


0.01205 


0.02716 


Cumulative 


0.00283 


0.00192 


0.00008 


0.00027 



TABLE VI. As in Table|V]but for background for pp -> e + e~p^ iss , with cuts for M H + = 100 GeV. 





VV Background 


tt Background 


Cuts 


Parton PYTHIA/PGS Parton 


PYTHIA/PGS 


Basic cuts 


0.42708 


0.33912 


0.58407 


0.40612 


Z pole veto 


0.74727 


0.73255 


0.86236 


0.85501 


150 GeV < M T2 


< 300 GeV 0.00260 


0.00196 


0.00000 


0.00000 


Cumulative 


0.00083 


0.00049 


0.00000 


0.00000 



TABLE VII. As in Table E] but for background for pp -)■ e + e"^ iss , with cuts for M H + = 300 GeV. 



M H+ 


- 100 GeV 


M H+ 


= 300 GeV 


Cuts Parton PYTHIA/PGS Parton PYTHIA/PGS 


Basic cuts 0.51713 


0.43680 


0.84810 


0.69845 


Z pole veto 0.90869 


0.90075 


0.97756 


0.97696 


Jet veto 0.68310 


0.57831 






M T2 cut 0.16875 


0.17320 


0.47802 


0.46847 


Cumulative 0.05417 


0.03941 


0.39632 


0.31966 



TABLE VIII. As in Table |V] but for the signal process pp —> fi + fi p™ lss via H + H . 



16 



VV Background tt Background 

Parton PYTHIA/PGS Parton PYTHIA/PGS 



Cuts 

Basic cuts 0.32959 
Z pole veto 0.73839 
Jet veto 0.62703 
M w+ < M T2 < 100 GeV 0.01324 
Cumulative 0.00202 



0.28226 0.52593 0.39048 

0.73098 0.86021 0.85489 

0.63204 0.01346 0.02282 

0.01554 0.00657 0.03675 

0.00203 0.00004 0.00028 



TABLE IX. As in Table [V] but for background for pp -> fi + ^-p^ iss , with cuts for M H+ = 100 GeV. 





VV Background 


tt Background 


Cuts 


Parton 


PYTHIA/PGS Parton 


PYTHIA/PGS 


Basic cuts 


0.32959 


0.28226 


0.52593 


0.39048 


Z pole veto 


0.73839 


0.73098 


0.86021 


0.85489 


150 GeV < M T2 


< 300 GeV 0.00288 


0.00239 


0.00000 


0.00000 


Cumulative 


0.00070 


0.00049 


0.00000 


0.00000 



TABLE X. As in Table|V]but for background for pp -)• fi + ^-p^ iss , with cuts for M H + = 300 GeV. 



larger than the QCD and parton density uncertainties on the VV and tt backgrounds; the 
overall cross sections of these backgrounds can also be normalized experimentally using Mt2 
regions below My/. 

For M H + = 100 GeV, the background after cuts depends sensitively on the shape of the 
background Mt2 distribution just above My/. This is controlled by the W width and the 
detector resolution for lepton momenta and missing p^', its shape should not suffer from 
QCD or parton-density uncertainties. For M H + = 300 GeV, the shape and normalization 
of the nonresonant tail of the VV background is especially important. This background is 
mostly Drell-Yan with an additional on-shell W boson radiated from one of the final-state 
leptons; the QCD corrections to such processes are well understood. Given enough statistics, 



M H+ 


= 100 GeV 


M H+ 


= 300 GeV 


Cuts Parton PYTHIA/PGS Parton PYTHIA/PGS 


Basic cuts 0.56131 


0.45743 


0.88249 


0.70832 


Jet veto 0.68783 


0.59528 






M T2 cut 0.16857 


0.16121 


0.47427 


0.46373 


Cumulative 0.06508 


0.04390 


0.41854 


0.32847 



TABLE XL As in Table W\ but for the signal process pp — > e ± /x =F p™ lss via H + H . 
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VV Background 


tt Background 


Cuts 


Parton 


PYTHIA/PGS Parton 


PYTHIA/PGS 


Basic cuts 


0.35835 


0.30423 


0.55297 


0.39556 


Jet veto 


0.65590 


0.68572 


0.01255 


0.02592 


M w + < M T2 


< 100 GeV 0.00860 


0.01207 


0.01585 


0.03018 


Cumulative 


0.00202 


0.00252 


0.00011 


0.00031 



TABLE XII. As in Table|V]but for background for pp ->■ e ± fi T p^ iss , with cuts for M H + = 100 GeV. 



VV Background tt Background 

Cuts Parton PYTHIA/PGS Parton PYTHIA/PGS 

Basic cuts 0.35835 0.30423 0.55297 0.39556 

150 GeV < M T2 < 300 GeV 0.00057 0.00049 0.00000 0.00000 
Cumulative 0.00021 0.00015 0.00000 0.00000 



TABLE XIII. As in Table|V]but for background for pp -> e ± fi T p^ iss , with cuts for M H + = 300 GeV. 

the shape of this background could also be normalized using the Mt2 region above M#+. 
Note also that the nonresonant tail of the VV background is significantly smaller for the 
e^fi^ final state than for the e + e~ and /i + yU~ final states, leading to a much higher signal 
purity in this final state for M#+ = 300 GeV as shown in the last line of Table IXIVI (for the 
lower charged Higgs mass this effect is swamped by the resonant-IT 7 contribution). 

The integrated luminosity required for a 5a discovery of H + H~ is displayed in Fig. [3] 
for M H + = 100 GeV and Fig. @] for M H + = 300 GeV, for each channel separately and 
for all three channels combined. We use the PYTHIA-PGS results and compute only the 
statistical significance. For the normal hierarchy with M#+ = 100 (300) GeV, we find 
5a discovery statistics with a minimum of 9 (56) fb -1 . For the inverted hierarchy, the 





Channel 


S/B Luminosity for 5a 






0.24 


78 fb- 1 


100 GeV 


/i-V~PT iss 


0.22 


88 fb" 1 






0.22 


40 fb" 1 




e + e-pf iss 


0.25 


526 fb" 1 


300 GeV 


fi+fi-pf™ 


0.25 


540 fb" 1 




e ± ^pf iss 


0.89 


73 fb" 1 



TABLE XIV. Signal over background (S/B) and luminosity required for a 5a discovery in a single 
channel for the three signal processes studied, for M H + = 100 and 300 GeV, assuming a degenerate 
neutrino spectrum so that BR(i7 + — > e + v) = BR(i/ + — > /x + i^) = 1/3. 
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FIG. 3. Luminosity required at the LHC (14 TeV) for a 5cr discovery if M H + = 100 GeV, for 
the normal hierarchy (NH, left) and inverted hierarchy (IH, right). The lines for each channel 
bound the range of required luminosities obtained by scanning over the 2a allowed ranges of the 
parameters of the neutrino mixing matrix and mass-squared differences. 




0.001 0.01 0.1 1 0.001 0.01 0.1 1 



lightest v mass [eV] lightest v mass [eV] 



FIG. 4. As in Fig. [3] but for M H + = 300 GeV. 



minimum is 8 (24) fb _1 . For the case of degenerate neutrino masses, 20 (57) fb^ 1 is needed. 
For degenerate neutrino masses, the luminosity needed for a 5cr discovery in each channel 
separately is given in Table IXIVl 
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IV. DISCUSSION AND CONCLUSIONS 



The two-Higgs-doublet model for Dirac neutrino masses studied here provides distinctive 
leptonic signatures at the LHC due to the characteristic decay pattern of the charged Higgs 
boson, controlled by the neutrino masses and mixing. We have shown that a simple set of 
cuts allows discovery of charged Higgs pairs with decays to i^p™ iSS with relatively modest 
integrated luminosity. In particular we found that a cut on the kinematic variable M?2 
provides very effective suppression of W pair and tt backgrounds for charged Higgs masses 
sufficiently above the W mass. 

In the inverted neutrino mass hierarchy, the large branching fractions of the charged 
Higgs to ev and \xv guarantees a 5a discovery for any allowed neutrino mass and mixing 
parameter values with only 20 (57) fb _1 for M#+ = 100 (300) GeV. The discovery potential 
remains remarkably good at = 300 GeV despite the rapidly falling charged Higgs pair 
production cross section because of the increasing separation of the signal M T2 distribution 
from the background. 

In the normal neutrino mass hierarchy, the large uncertainty on the neutrino mixing angle 
823 leads to parameter regions in which the charged Higgs decays predominantly to tu, with 
a branching fraction to light leptons below 40%, resulting in poor discovery sensitivity in 
the light lepton channels studied in this paper. Away from these parameter regions, the 
discovery prospects are only slightly worse than in the inverted hierarchy. 

As more stringent experimental limits are placed on the neutrino parameters from neu- 
trino oscillation experiments and direct searches for the kinematic neutrino mass in beta 
decay, the predictions for the charged Higgs branching ratios in this model will tighten. 
For example, one goal of the currently-running T2K long-baseline neutrino oscillation ex- 
periment in Japan is to improve the measurement accuracy of sin 2 (2^23) by an order of 
magnitude [2l|, which would reduce the 2a spread in the charged Higgs branching ratios to 
\w and rv at low lightest-neutrino mass from the current ±30% to about ±10%. Sensitivity 
to the neutrino mass hierarchy relies on detection of a nonzero #13, a major goal of T2K and 



22|. The 



the longer-baseline U.S.-based experiment NOi^A currently under construction 
ratios of the signal rates in the three channels considered here would allow the normal, 
inverted, and degenerate neutrino spectra to be differentiated, providing a key test of the 
connection of the model to the neutrino sector. 
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Measurement of the charged Higgs branching fractions will also provide some sensitivity 
to the mass of the lightest neutrino. For a lightest neutrino mass between about 0.01 
and 0.1 eV, the charged Higgs branching ratios vary dramatically with the value of the 
lightest neutrino mass (Fig. [T]); once the measurement of #23 from neutrino oscillations has 
improved, measurement of the ratio of the ev and \xv modes will provide sensitivity to 
the lightest neutrino mass in this range. This is nicely complementary to the prospects 
for direct kinematic neutrino mass determination from the tritium beta decay experiment 



KATRIN, which is designed to be sensitive down to neutrino masses of about 0.2 eV [23 1 — 
i.e., at the lower end of the degenerate part of the spectrum — and is scheduled to begin 
commissioning in 2012 [24]. We note that, because the neutrinos in this model are Dirac 
particles, neutrinoless double beta decay experiments will have no signal and will thus not 
be sensitive to the neutrino mass scale. 

The mass of the charged Higgs is also accessible at the LHC through the signal event 
kinematics. In particular, the signal Mj> 2 distribution is flat up to an endpoint at the charged 
Higgs mass, as shown in Fig. [2J A fit to this distribution on top of the background should 
provide a measurement of the charged Higgs mass. This would allow a valuable cross-check of 
the charged Higgs pair production cross section together with the visible branching fractions 
as predicted by the neutrino parameters. The pair production cross section is sensitive to the 
isospin of the charged Higgs through its coupling to the Z boson, allowing the two-doublet 



nature of the model to be established [5]. 

We finally comment on the applicability of our results to two other neutrino mass models 
that contain a charged Higgs boson. First, the Z 2 model of Ref. j^j contains a charged Higgs 
with partial widths to leptons and LHC production cross section identical to those in our 
model. The charged Higgs in the Z 2 model differs from ours in that it can also decay to 
W + a, where the neutral scalar a is extremely light due to the spontaneous breaking of the 
Z 2 symmetry. This competing mode dominates unless H + is not much heavier than Myy and 
the neutrino Yukawa couplings are 0(1) (4] (this parameter region is forbidden by standard 
big-bang nucleosynthesis, but the Z 2 model already requires nonstandard cosmology due 
to the very light scalar a). For this parameter range, then, our results should carry over 
directly. For smaller Yukawa couplings or a heavier H + , the decays to leptons used in our 
analysis are suppressed, resulting in a smaller signal on top of the same background. 

Second, neutrino masses of Majorana type can be generated by the so-called Type II 
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seesaw mechanism 



25| . in which an SU(2)-triplet Higgs field X = (x ++ ,x + ,x°) with 



very small vev is coupled to a pair of SM lepton doublets. LHC phenomenology for this 
Higgs-triplet model was studied in Ref. 26), which considered signatures from x ++ X an d 
X ++ X~ (and the conjugate process) at the LHC. While the decay branching fractions of 
X + in this model are identical to those of the charged Higgs in our Higgs-doublet model, 
the LHC production cross section for X^X~ m the triplet model is about 2.7 times smaller 
than for H + H~ in the doublet model [5|, due to the different isospin of x + which modifies 
its coupling to the Z boson. The signals studied here would thus have a S/B of less than 
10% for most channels, potentially leading to problems with background systematics. For 
sufficiently high charged Higgs mass, though, the /i ± e T channel would still have a decent 
S/B (33% for M H + = 300 GeV and a degenerate neutrino spectrum); the reduced cross 
section in the triplet model would however require an integrated luminosity close to 300 
fb _1 for discovery. In any case, searches for the doubly-charged scalar would yield an earlier 
discovery of the triplet model. 
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